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ABSTRACT 

We propose that the intensity changes and spectral evolution along the M87 
jet can be explained by adiabatic changes to the particle momentum distribution 
function and the magnetic field. This is supported by the lack of any significant 
variation in the radio-to-optical spectral index along the jet and the moderate 
changes in radio brightness. Assuming a simple scaling law between magnetic 
field and density, we use the deprojection of a 2 cm VLA intensity map by 
Sparks, Biretta, & Macchetto (1996) to predict the spectral evolution along the 
jet. 

We derive limits for the magnetic field and the total pressure by comparing 
our results with the spatially resolved fit to spectral data by Neumann, 
Meisenheimer, & Roser (1997) of a model spectrum that cuts off at ^ lO^^Hz. 
To explain the weakness of synchrotron cooling along the jet , the magnetic field 
strength must lie below the equipartition value. Although the inferred pressure 
in the limit of nonrelativistic bulk fiow lies far above the estimated pressure 
of the interstellar matter in the center of M87, bulk Lorentz factors Fjct in the 
range of 3 — 5 and inclination angles ^los ^ 25° lead to pressure estimates close 
to the ISM pressure. The average best fit magnetic fields we derive fall in the 
range of 20 — 40 /xG, departing from equipartition by a factor 1.5 — 5. 

This model is consistent with the proposal by Bicknell & Begelman (1996) 
that the knots in the M87 jet are weak, oblique shocks. First-order Fermi 
acceleration will then have a minimal effect on the slope of the radio-to-optical 
spectrum while possibly accounting for the X-ray spectrum. 
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1. INTRODUCTION 

The jet in M87 has been observed in many wavelength bands, yet some puzzles remain 
about the nature of its emissivity. Polarization observations in both the optical and radio 
have shown that the emission from the jet at wavelengths longward of ~ 100 A is most 
likely of synchrotron origin (|Baade 1956| , |Owen et al. 1989|) . Collimated structure can be 



traced back to O.Olpc from the core using VLBl ( |Junor fc Biretta 1995|) , and terminates 
25" away from the core (2 kpc at the assumed distance of 17 Mpc) in the western radio 
lobe, with optical emission still detectable at this distance. The jet looks very similar at 
optical and radio wavelengths, although Sparks, Biretta, & Macchetto (1996, hereafter 
showed that differences (e.g., in the transverse brightness profiles) do exist. 



Far from having a smooth appearance, the jet exhibits a series of bright knots at 
intervals roughly 2". 5 apart. The nature of the knots is uncertain, but they are usually 
attributed to internal shocks from either flow instabilities (Bicknell & Begelman 1996, 
hereafter pB| ) or variable outflow at the source ( [Rees 19"78D . Both the knots and the 



interknot regions exhibit power law spectra of index ~ 0.5 in the radio, connecting to 
the optical data with a power law of index a-^o ~ 0.65, and steepening to ao between 1.2 
and 1.8 in the optical. This steepening trend is also found in observations at infrared (e.g., 
Stocke, Kieke, &: Lebofsky 1981| ) and ultraviolet ( [Ferola fc 'I'arenghi 1980| ) wavelengths, as 



confirmed recently by i/ST observations ( |Boksenberg et al. 1992| , [SB1V1| ). X-ray observations 



made with Einstein Observatory ( [Biretta, Stern, fc Harris 199 1| ) and ROSAT ( [Neumann et 



|al. 1997a| and [Reynolds et al. 1996b| ) reveal X-ray emission from several spots along the jet 



(mainly the core and knot A, possibly also from knots D and B). However, the origin of the 
X-ray emission is unknown and it is not clear whether the spectrum breaks between optical 
and X-ray wavelengths to a spectral index of aox ~ 1-4, with the X-ray emission still being 
of synchrotron origin, or whether the X-ray emission is produced by a different mechanism. 
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e.g., inverse Compton scattering or bremsstrahlung. 

The onesidedness of the M87 jet has been interpreted as the result of relativistic 
beaming of the emission from an intrinsically bipolar jet, which implies bulk velocities 
corresponding to Lorentz factors of Fjct ^ 2 with line of sight inclinations of 6'los ^ 30°. 
This would explain the existence of the radio lobes on both sides of the core. New 
infrared detections of counterjet emission (Stiavelli, Peletier, & Carollo 1997) seem to 
confirm this hypothesis. The Doppler beaming interpretation is bolstered by proper motion 
measurements of the knots ([Biretta, Zhou, fc Owen 1995|) , which show characteristic 



velocities fproper ~ 0.5 c in knots A and B, with some features exhibiting much larger proper 
velocities (a subfeature of knot D appears to show superluminal motion). These motions 
broadly support the interpretation of the knots as relatively weak, oblique shocks moving 
down the jet with pattern speeds significantly smaller than the bulk speed of the fiow 



One might hope to detect the effects of synchrotron cooling and relativistic particle 
acceleration by studying the spectrum as a function of position along the jet. Such 



measurements (^BM|) show that the radio-to-optical spectral index, a^o, is very nearly 
uniform, while the optical spectral index, ao, is anti-correlated with the brightness, i.e., the 
optical spectrum is flatter in regions of higher intensity. Modeling the optical steepening 
as a high-energy cutoff imposed on a power law spectrum, Neumann, Meisenheimer, & 
Roser (1997, hereafter [NMRJ) and [Meisenheimer, Roser, fc Schlotelburg (1996)| flnd a 



corresponding correlation between brightness and cutoff frequency, i.e., a higher frequency 
at higher intensities. Both these results are striking in the lack of a strong secular decline 
in the cutoff frequency with distance from the core, as would be expected naively if the 
steepening were due to synchrotron coohng. 

Indeed, these observational results do not compare well with simple quantitative 
models of synchrotron cooling in the M87 jet. The usual assumption of equipartition leads 
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to estimates of the magnetic field of order 300 fiG, with values up to 500 /xG in knot A. 
The synchrotron lifetime for electrons with Lorentz factors of 7 ~ 10^, needed to produce 
the optical emission in a 300 fiG field, is only 2.3 x 10^'^ sec, which, for mildly relativistic 
bulk velocities (e.g., 0.5 c) implies a travel distance of less than 120 pc. But the projected 
length of the jet is about 2 kpc, and even the distances between the most prominent knots 
are longer than the estimated cooling length. Yet the spectrum between radio and optical 
bands remains remarkably constant along the jet, with only minor variations in the optical 
spectral index. 

This presents a paradox: After a travel distance of 2 kpc one would expect the cutoff 
to have moved to below 10^^ Hz, a factor of 10'^ smaller than what is actually observed. The 
discrepancy becomes worse if the magnetic field is stronger than the equipartition value, 
as suggested by Owen, Hardee, & Cornwell (1989) in order to explain the confinement of 
the overpressured jet via the magnetic tension force. (An alternate explanation for the 
confinement — that the radio cocoon surrounding the jet is overpressured with respect to 
its surroundings — has been proposed by PB| .) Three explanations for the discrepancy have 
appeared in the literature: 

1. First-order Fermi re-acceleration in the knots, interpreted as shocks, could produce 
high energy electrons (and possibly positrons) from a synchrotron-cooled distribution, 
with a power law of roughly the index observed. It could also explain the observed 
X-ray emission. However, since the power-law index produced by Fermi acceleration is 
a strong function of the compression ratio one would need some fine tuning to explain 
the observed constancy of the radio-to-optical spectral index, which does not equal 
the limiting value for a strong adiabatic shock. Furthermore, there does not seem 
to be a significant amount of cooling between the knots even at the highest optical 
frequencies observed, which would be expected for the assumed S-fields and interknot 
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distances. 

2. Particles could be transported in a loss-free channel in the interior of the jet, with the 
bulk of the emission produced in a thin outer layer of high magnetic field strength 
( Pwen et al. 198"9| ). In this picture, the knots and filaments would be interpreted as 
instabilities with greatly increased magnetic field strength wrapped around the jet. 
The emission would then be fed by particles from this channel. The radio brightness 
profiles across the jet seem to suggest a limb-brightened emission, but a reinvestigation 
of the HST observations ( ^BM] ) shows that the optical emission is more concentrated 



to the inner regions of the jet. Also, |SBM| 's deconvolution places the brightest spots 



in the jet interior. This argues against a field-free zone in the jet interior. 

3. On-the-spot reacceleration by a yet unknown process could maintain the cut-off 
particle momentum at the observed level, as has been proposed by Meisenheimer et 
al. (1996) in a model similar to ours (see § [4.2|) . This has the advantage of explaining 
all the observed features, but invokes unknown physics to explain the apparent lack 
of cooling. 

Inspired by the observed correlation between the emissivity variations and the cut-off 
frequency, by the newly deconvolved volume emissivity (|SBM| ), and by the new evidence for 
relativistic bulk velocities (|Stiavelli et al. 1997|) , we propose a simple way of explaining the 



observations. The only standard assumption we give up is the assumption of equipartition, 
which does not seem to have a very firm physical foundation anyway. Magnetic fields 
smaller than equipartition by a factor of 1.5 — 3, coupled with bulk Lorentz factors in 
excess of 2 — 3, can readily explain the lack of evidence for synchrotron cooling. In our 
model the fluctuations of the cutoff frequency are produced by weak shocks, so that the 
influence of the compressions on the plasma distribution function can be considered to be 
adiabatic. As a result, we are able to explain the general behavior of the cutoff reasonably 
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well; the magnetic fields we derive from our own fits to the data are below equipartition 
with a convincing level of confidence, but the inferred total pressures do not necessarily 
need to exceed the equipartition values. Additionally, for relativistic jets the equipartition 
values for i?-field and pressure are less than in the nonrelativistic limit, so the pressures we 
derive can fall below the equipartition value in the nonrelativistic case. 

We organize the paper as follows. In § |^ we present a nonrelativistic treatment of the 
synchrotron emissivity, taking into account cooling but assuming that the particles and 
fields respond adiabatically to changes in the flow density. Once relativistic effects are 
incorporated into the treatment in § ^, we use the data of |SB1VI| and |M 1V1K| to constrain the 



magnetic field strength of the jet. Section ^ discusses confinement and stability of the jet in 
the light of the pressures derived from the results of § |^, the production of X-ray emission 
in knot A, polarization, and limits on the particle acceleration site; and § |^ gives a brief 
summary of the results and future prospects. 



2. ADIABATIC EFFECTS ON SYNCHROTRON EMISSION 

Our model rests on the hypothesis that Fermi acceleration is unnecessary to explain the 
fluctuations of radio-to-optical emissivity and cutoff frequency along the M87 jet. Given 
certain assumptions about the orientation and degree of disorder in the magnetic field, and 
the degree of anisotropy permitted in the relativistic electron distribution, we can relate 
changes in both the emissivity and the cutoff frequency uniquely to changes in the density 
of the jet fluid. These adiabatic effects are readily combined with the effects of synchrotron 
cooling ([Coleman fc Bicknell 198"8|) . In effect, given the emissivity map of SBM| , we can 



predict the run of the cutoff frequency along the jet, and vice-versa. Since we also have 
NMR's observations of the cutoff frequency as a function of position, our adiabatic model 
is subject to a powerful self-consistency check. 



- 8 - 



Observationally, the main changes in emissivity and cutoff frequency are rather 
locahzed, and associated with the positions of the knots. These small-scale fluctuations 
therefore provide the strongest check on our assumption of adiabaticity. The large-scale 
trends then determine the best fit to the magnetic field strength. We have already seen 
that the apparent lack of a large-scale synchrotron cooling trend is incompatible with a 
magnetic field strength as large as the mean equipartition value (see Fig. |T]). Neglecting 
relativistic and projection effects, it would require a magnetic field as low as 25 to 
obtain a cooling length of 2 kpc at optical frequencies. This is an order of magnitude 
smaller than the mean equipartition field and would require a total (particle + magnetic) 
pressure of 1.2 x 10^^ dyn cm^^ to produce the observed average amount of synchrotron 
emission, compared to an equipartition value of ~ 4.0 x 10~^ dyn cm~^. 

As noted earlier, the physical basis of equipartition is weak. Estimates of -B-fields 
in radio hot spots based on the synchrotron cooling time indicate that in some cases 
equipartition might be correct up to a factor of ~ 2 ( pVIeisenheimer et al. 1989| ), but the 



conditions in the jet might very well be different from those in the lobes. We are therefore 
free to consider the magnetic field strength to be a free parameter. Applying our adiabatic 
model to the observational data, we can derive an estimate for the magnetic field strength. 
In the case of the M87 jet, this estimate lies below equipartition, even when relativistic 
effects are taken into account (§ |]). 

Our model for the evolution of the particle distribution function follows that of Coleman 
& Bicknell (1988). We assume that pitch angle scattering due to plasma micro-instabilities 
keeps the particles close to an isotropic distribution in the fluid rest frame. In the absence 
of cooling, this would imply that the relativistic electrons respond to compressions like a 
7adiabatic = 4/3 (i.e., ultrarclativistic) fluid, but this behavior will be modified by synchrotron 
cooling. Because the magnetic field is frozen into the plasma, its strength should change as 
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the plasma density fluctuates along the jet. 

Depending on the orientation and the degree of disorder of the field, its variation will 
depend roughly on the density change to some power (: B{r) oc g{r)'^, where g is the proper 
particle density and r is the distance from the core. For a completely disordered magnetic 
field, C = |) whereas for a homogeneous field the power depends on the orientation of the 
field with respect to the compression normal, with ( = 1 for an orthogonal orientation, 
C = for parallel orientation. Since the polarization of the jet is of order 10%-20% in the 
interknot regions, compared to the maximum polarization of 70% for a homogeneous field, it 
is likely that the magnetic field has a disordered component, so it is reasonable to assume ( 
to be of order | (but see §4.4] ). The ordered 5-component is aligned with the jet axis almost 
everywhere except in the brightest knots, as can be seen from polarization measurements 
( Piretta fc Meisenheimer 1993| ). We do not know its orientation with respect to the 



compression normal, because the orientations of the (presumably oblique: see pB|) shocks 
are unknown. However, our results are not very sensitive to what the actual value of ( is, 
as we will show later. We therefore make the simplifying assumption of a single exponent 
describing the field variations. Using the scaling relation for the synchrotron emissivity of 
a power law momentum distribution (/(p) d^p oc d^p oc dp, corresponding to a 



spectral index of a = ^y^) under adiabatic compression (e.g., [Coleman fc Bicknell 1988| ): 

J oc (5sin^9)^+"^^i/-", (1) 

where is the angle between magnetic field and line of sight, we can express the field 
relative to its value at tq = 0".5, the (arbitrary) injection point at which we start the 
calculation, as a function of the emissivity ratio j/j(ro): 

B = B{ro)ij/jiro)f (2) 

where ^ = [1 + a + (2a + 3)/(3C)]" . A necessary condition for this approach to be valid 
(in addition to the assumed isotropy and the absence of Fermi acceleration) is the assumed 
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steady-state injection of relativistic particles and fields by the central engine, which allows 
us to relate densities and fields at each r to the corresponding values at the injection point 
vq at a given instant of time. 

Is it plausible to neglect Fermi acceleration in the shocks that comprise the knots? 
Except for knot A, the brightness changes along the jet are moderate. In knot A, the 
brightest feature, the emissivity changes by a factor of order 10 ( pBM|) , which, if entirely 
due to a sudden compression of the plasma, can be produced by a proper compression ratio 
of r ^ 2.7 as measured in the respective rest frames of the plasma if we take C = §• For the 
other knots we infer smaller density contrasts, r ^ 1.5. Consistent with this observation, 
we will henceforth take the knots to be weak shocks, in accordance with the suggestion by 



BB| that the knots are highly oblique (and therefore weak) shocks. Thus, because Fermi 
acceleration leaves the spectral index unchanged if the shock is weak enough, we will 
henceforth neglect its effect on the cutoff frequency. Section |^ discusses Fermi acceleration 
in more detail, with particular focus on the possibility of Fermi acceleration occurring in 
knot A. 

Furthermore, because the shocks are believed to be oblique, we take the fluid velocity to 
be constant to first order, both in magnitude and direction. For the shock jump conditions 
in the non-relativistic limit (which we consider in this section) the velocity component 
perpendicular to the shock plane v± is inversely proportional to the density, thus for a 
proper compression ratio of 2.7, as seen in knot A, the perpendicular velocity component 
should change by a factor of 0.37. For highly oblique shocks, v± is small compared to f y 
and the velocity will not change significantly. Moderate changes in velocity would be easy 
to incorporate in principle; yet with our current ignorance of the velocity field and shock 
parameters, such a level of detail is unwarranted. We will comment on the validity of this 
assumption in § pi. We also postpone a treatment of the motion of the knots until § El, and 
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assume them to be stationary for the rest of this section. 

With these assumptions we are ready to calculate the downstream cutoff frequency for 
a given initial cutoff momentum in the injected particle distribution and a given i?-field at 
Tq. We use the transport equation as presented by Coleman & Bicknell (1988): 

written in the rest frame of the fluid. Here, f = f{p) is the electron distribution function 
and A = ^^^B'^ is the synchrotron loss term. The equation is valid for the assumed case of 
isotropy and negligible Compton losses (for a brief discussion of Compton losses see 
The solution of this equation is 

f{p{r)) = Mpo) (^] (4) 

where /o is the injected momentum distribution and 

(g/go)^Po 

:i+poi:i:!)Ait')ig{t')/gor^df) 

(see [Coleman fc Bicknell 1988|) where the subscript denotes the values at injection point 



?('■) = - — (s) 



Equation (^ describes how the momentum of a given particle changes along a 
streamline. Thus, if the distribution initially cuts off at pc^, we can calculate the cutoff 
momentum Pc{r) downstream. Because in our model the density g{r) is proportional to 
B(r)( and because we know the scaling of B with r from equation (^, we can eliminate g 
and B/Bq from equation (^). The remaining parameters are Bq, p^fi, and j(r)/jo, the latter 
being provided by the Sparks et al. data. 

The cutoff momentum p^ is related to the observed cutoff frequency by the expression 

3e 
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Equation (^) contains another parameter, ■(9, the angle between the hne of sight and the 
magnetic field. For now we shall set the factor sin^ ^ 1, which is valid for disordered fields, 
since 1) the regions in which the field is perpendicular to the line of sight have the highest 
emissivity, and 2) assuming randomly oriented fields, half of the field orientations lie in the 
range from 60° to 90° to the line of sight, i.e., sin?? > 0.866. Thus the cutoff frequency is 
mainly determined by field orientations close to 90° or sin^ ^ 1. 

We can now determine the free parameters Bq and Pc,o by applying a least chi-squared 
method to fit the observed cutoff frequency z/c,obs with the value determined from equations 
(|) and (^, using the emissivity map j{r)/jo provided by |SB1V1[ We prefer to average 



the emissivity across the jet, which minimizes small scale variations probably due to the 
deprojection procedure. Because [NMR| also averaged across the jet, this seems to be the 



most appropriate way of calculating the cutoff frequency. Figure (calculated for a bulk 
Lorentz factor of Fjet = 1.1, a radio spectral index of c^r = 0.5, 6'los = 90°, and C = |) 
shows the observed cutoff frequency z/c,obs (vertical bars) with error bars and the best fit 
curve (solid line), which seems to reproduce the scaling of the cutoff frequency reasonably 
well. (The radio spectral index seems to break to a ~ 0.65 at ~ 10 GHz, so we have used 
both a = 0.5 and 0.65 in our fits with insignificant differences in the average parameters 
but smaller chi-squared for a ~ 0.65; see §^) For comparison the plot also shows the best 
fit cutoff frequency for equipartition i?-fields (dashed line). The mean 5-field is of order 
10 fiG, even smaller than the zeroth order estimate made at the beginning of this section. 
Assuming (arbitrarily) a lower cutoff at z/ = 10'' Hz and the observed high-frequency 
cutoff at z/ ^ 10^^ Hz yields an average total pressure of 8 x 10^^ dyn cm^^ for the given 
parameters, compared to an equipartition value of p ~ 3 x 10^^ dyn cm~^. In calculating 
absolute values for both pressure and 5-field, projection (i.e., foreshortening and length 
scale) effects must be taken into account, since the emissivity was derived for a side-on 
view of the jet. This introduces a factor of sin 6'los in intrinsic emissivity and pressure for a 
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given magnetic field and a factor of (sin6'Los)^ in equipartition pressure. 

It seems that the proposed modest compressions can account for the fiuctuations seen 
in the spectral cutoff, and the large scale decrease in Uc is well reproduced. It would be 
helpful to determine both the emissivity and the spectral index maps from the same method 
and data, thus eliminating errors due to different reduction procedures. We comment on 
the deviations and uncertainties in this fit in § ^. It is important to note that this technique 
should be independent of what the actual shape of the particle distribution is, because it 
simply tracks the behavior of a single feature in the spectrum, which could be identified 
with either a break or a cutoff. 

As seen in this section, one runs into problems with the jet pressure for nonrelativistic 
bulk velocities. Also, the observed mildly relativistic proper motion of the knots and the 
jet's onesidedness favor a relativistic interpretation, as do the knot spacing and morphology 
( [BBp . The results indicate that Lorentz factors of order 2 — 5 fit the observations best. In 
the next section we will investigate the effects of these suggested relativistic bulk velocities. 

3. RELATIVISTIC EFFECTS 

Relativistic motions not only explain the onesidedness of the M87 jet, but also help to 
solve the synchrotron cooling problem mentioned in the introduction. The travel time in the 
electron rest frame is reduced by a factor Fjet due to time dilation, the intrinsic emissivity 
is reduced by a factor (where D is the Doppler factor D = [rjet(l — [3 ■ cos6'los)]~^)5 

and the intrinsic cutoff frequency is Doppler shifted downward by a factor D. As a result, 
the apparent synchrotron lifetime can be a significant underestimate of the intrinsic value. 

Biretta (1993) estimates the lower limit on the jet-to-counterjet radio brightness ratio 
to be > 150 — 380 (the higher value corresponds to the assumption that jet and counterjet 
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have identical appearance). Based on this limit we adopt line-of-sight angles < 35° and 
Lorentz factors > 2. Note that even Doppler factors smaller than unity can lead to a large 
jet-to-counterjet brightness ratio, as the counterjet brightness is severely reduced for large 

Tjet- 



We repeat the analysis of § ^ using the emissivity profile of |SBM| , this time corrected 



for Doppler boosting and projection effects. Again, for fitting the cutoff frequency only 
emissivity ratios are important, so these corrections do not change the fitting procedure 
as long as changes in the bulk velocity can be neglected. We are confident that at least 
the direction of the flow is not changed significantly before knot C. The only knot for 
which such effects could be important is knot A, because it displays a jump in emissivity 
of 11, whereas in the other knots the emissivity is increased by a factor of order 3 only, 
which implies very moderate compressions. Using the relativistic continuity equation for an 
oblique shock we have estimated the post-knot A Lorentz factor to be rjet,A+ ~ 3 for a 
pre-knot A rjet,A- = 5, a compression ratio of 3 and intrinsic obliquities ~ 60°. Although 
the change in Fjei might seem large at first glance, the impact such a velocity change has 
on our fits is not large, as we will explain below. We modify equation (|^) to follow the 
electron distribution in the fluid rest frame by replacing t with r and with where r 
is the proper time. The same changes apply to equation (^. Strictly speaking, the fluid 
frame is not an inertial frame and we would have to include accelerational terms into the 
equation, introducing an anisotropy. But our assumption should be adequate, provided that 
isotropization takes place over short enough scales. 

Treating the response of the distribution function to compressions as adiabatic and 
assuming isotropy (i.e., an adiabatic index of |) we calculate the changes in B and cutoff 
momentum (measured in the fluid frame) from the emissivity changes in the fluid frame. 
The substitution t —* t takes care of the time dilation effects. 
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In order to incorporate the observed motion of the knots, we must correct A{t') and 
g{T') in equation (j^) for hght travel time effects between the source and the observer. This 
is because the knots move during the time it takes a particle to travel from tq to r. We need 
to know the ratios B / Bq and q/ Qo experienced by a particle as a function of proper time 
t' in order to be able to do the integration in equation (|^). Because we infer the relative 
values of B and ^ at a given position from the emissivity ratio j(?"fluid)/jO) it is important 
to know the velocity of the emissivity pattern. 

We assume that the pattern of density and field fluctuations retains its shape and 
moves along the jet at a fixed speed ^pattern, taken to be smaller than t>fluid and set equal to 
0.55 c everywhere in our calculations for simplicity. If the present field distribution (i.e., at 
a given time t = in our frame, corrected for light travel time effects) is expressed by B{r), 
then the field at time —t' and position r' is given by B{r' + fpattom^')- Now, for a particle 
currently at r, the equation of motion is r' = r — ffluidi^'- Therefore, the field distribution 
experienced by the particle as a function of time is i?[r — (f fluid — ^^pattern)^']- Appropriate 
modifications to equation @ are straightforward. The effect of the pattern speed on the 
result is not very dramatic, reducing xLin about 6%. 

With this set of assumptions we can once again proceed to integrate the modified 
equation (|^) for various 6'los and Fjet- Using a minimum chi-square routine we can 
determine the best-fit values for Bq, and po. Relation (||) then yields B{r). 

Figure ^ shows a chi-square plot for Fjct = 3, 6'los = 25°, and C = |- The equipartition 
value for the average 5-field is shown as a shaded area at 89 fiG. The upper limit on -Bmcaru 
set by 2x^jjj contours, lies at 49 fiG, 75% above the best fit value, -Bmcan = 28 /iG. The 
lower limit set by 2%^;^ is ~ 5 fiG, 80% below the best fit value. The average equipartition 
field of i?moan ~ 89 fiG lies above even the 5Xmin contour. The lower limit on B is not 
nearly as strict, due to the fact that cooling is not dominant, i.e., we can produce a similar 
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spectral behavior by reducing the magnetic field and increasing the particle energy, which 
produces the tear-shaped appeareance of the contours. Thus, strictly speaking, the best fit 
values for the i?-field should be regarded as upper limits. 

The reduced Xmin values (i.e., Xmin divided by the number of degrees of freedom) fall 
above 44, which is uncomfortably high. However, because we do not have formal errors 
for the emissivity deprojection by |SBM| , which will introduce a significant uncertainty, a 



high value for Xmm that discouraging. Estimating the average uncertainty in 

the emissivity by comparing the averaged emissivity to that derived from taking only a 
slice along the jet yields an uncertainty of order 50%, which leads to uncertainties in the 
predicted cutoff frequency of roughly 20%. This is significantly higher than the formal error 
in [NMR] 's data and will reduce the Xmin ^ factor of approximately 10. 

The Xmin values are dominated by the region beyond knot A. The post-knot A residuals 
in our fit are not larger than the residuals in the pre-knot-A region, but because the post 
knot A region is brighter, the error bars on the measured cutoff frequency are smaller, 
which increases the Xmin- The deprojection procedure, which assumed an axially symmetric 
fiow, breaks down beyond knot A, which will introduce significant uncertainty. Also, 
non-uniformities in the emissivity could lead to large errors if the optical emission peaks at 
different locations than does the radio emission. Field orientation effects and changes in Fjet 
and 6'los might also contribute to the error. We performed the same procedure just out to 
knot A and found that, with the same parameters, the reduced Xmin shrinks to 13. Leaving 
B and Pc as free parameters reduces Xmin 10, but also reduces the i?-field significantly. 
Because in this case the algorithm mainly fits the region around knot A (where the error 
bars are smallest), we cannot expect the global run of Uc to have significant impact on 
the fit, which would be necessary to extract information about the average magnetic field. 
We conclude that the reproduction of fine detail is not satisfactory in the region beyond 
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knot A. However, the gross run of z/^, which is principally responsible for constraining our 
parameters, is reasonably well reproduced. 

The best-fit average magnetic field -Bmean, as a function of Fjet and 6'los5 is plotted in 
Figure |^. Figure || shows the average ratio (B/Bcq) as a function of Fjet for ^^los = 15° 
to 30° in increments of 5°, and the area corresponding to the limiting jet-to-counterjet 
brightness ratio of 150 - 380. Note that the equipartition magnetic field has to be corrected 

2 + a 1 

by a factor of D 3+q (^sin 6'los) for projection and Doppler boosting of the emissivity; 
this has already been taken into account in the figure. Clearly, for Fjet in the range 3 — 5 
and 6'los ^ 25° the departure from equipartition is not very large (roughly a factor of 



is small compared to the expected errors introduced by the simplifications we made and 
to the range in B allowed by our minimum chi-square procedure, at most 12% for C2 = 
and small Fjet. This is not a very reasonable value for ( in any case, because the field has a 
random component, thus ( should be higher, and the probability of the field being in the 
shock plane (thus having ^ = 1) is twice as high as for the field being normal to the shock. 
We conclude that our ignorance of the precise behavior of B under compression is not a 
serious obstacle to the application of our model. 

We also tested the impact a change in Tjet at knot A might have on our results. 
As we mentioned earlier, the best fit i?-field values we derive are upper limits. This 
is the reason why a change in Tjet at knot A does not change our results significantly: 
generally, lower Lorentz factors require lower fields to explain the observed lack of cooling. 
If the jet is slowed down beyond knot A, we will need lower average fields to fit this 



0.2 < (fi/ficq) < 0.6). 



In order to test the dependence of the best fit Bq on the parameter ( we have 
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region. However, lowering the field does not change the quality of the fit much (the is 
essentially unchanged), so the global field strength is simply set by the region with the lower 
Lorentz factor, rjet,A+ (the relative scaling of B is still determined from equation 0, taking 
relativistic beaming into account). We have introduced by hand a change of Fjei at knot 
A into our model (we solved the continuity equation at the shock, assuming an obliquity 
of 60°, for the velocity change that would reproduce the observed emissivity jump of 11, 
including relativistic beaming and adiabatic compression), and calculated the fractional 
deviation ( ^ | = l-^""'fo™~-^brcak ^j-^g derived averaged i?-field. Here, -Buniform is the best 
fit average i?-field derived for uniform Tjet and -Bbreak is the best fit field for a jet slowing 
down from Tjet,A- to Tjet,A+ at knot A. Figure |^ shows (■^)p for uniform jet models with 
Tjet set to either Tjet,A- or Tjet,A+ (filled light and dark grey regions, respectively). The 
latter is always less than 18% for the parameter range we used. Note that for post-knot-A 
Tjet.A+'s above 5, the pre-knot-A Tjet,A- exceeds 8.5, thus Tjet,A+ ^ 7 can be ruled out on 
the basis of gross energy balance arguments (see next section). 

For completeness we have shown the deviation of the best fit average 5-field 
f-^) = l-^"! -^"2 1 fQj, 2 cm radio spectral index of a2 = 0.65 instead of ai = 0.5 as the 
black region in Figure ^ One can see that the difference is negligible compared to other 
uncertainties. 



4. DISCUSSION 

In the preceding sections we demonstrated that a) magnetic fields slightly below 
equipartition and b) moderately relativistic effects are able to explain the general behavior 
of the spectrum in M87. In this section we will examine the confinement properties of the 
jet and compare our model with a previous model by Meisenheimer et al. (1996). We also 
comment on the production of X-ray emission in knot A, and on the consistency of our 
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model with polarization measurements, and we estimate the minimum distance from the 
core at which particle acceleration has to occur. 



4.1. Confinement 

Naturally the question arises whether the jet can be confined under the conditions we 
proposed above. The usual assumption for a jet to be confined is that it is in pressure 
equilibrium with its surroundings. Alternatively, one could imagine the jet to be freely 
expanding into an underpressured surrounding medium. 



BB argue that in order to produce shocks via Kelvin-Helmholtz instability, some 
interaction between jet and surrounding medium has to take place, as opposed to a free 
expansion scenario. They also show that the minimum Lorentz factor Fjet for a freely 
expanding jet with no cold matter content is at least 13, much higher than the values we 
have used above. Because we would most certainly fall out of the beaming cone for such a 
high Fjet jet, the intrinsic emissivity would be much higher than the observed value. As 



BB] point out, the energy fiux of the jet would far exceed the estimates made on the basis 
of the expanding bubble the M87 jet blows into the ISM. A Fjet that high would also raise 
questions about the location at which the jet is decelerated to nonrelativistic velocities, and 
seems inconsistent with the claimed detection of IR counterjet emission by Stiavelli et al. 
(1997). We can therefore rule out the picture of a freely expanding jet. As a consequence 
we need a mechanism to provide confinement, i.e., we need to set the jet pressure in relation 
to the ambient pressure. 



The ambient gas pressure in the center of M87 has been derived by [White fc Saraziii 



1(1988)1 from fitting cooling fiow models to the Einstein X-ray observations. The values 
they find fall into the range pisM = 1 x 10"^'' dyn cm~^ to pisM = 4 x 10~^° dyn cm~^. 
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It is important to note that the pressure of the interstellar medium in M87 might not 

's 



be representative of the pressure of the immediate environment of the jet. In fact, [BBf 
analysis of the helical Kelvin-Helmholtz instability leads to the conclusion that the ambient 
medium of the jet is significantly overpressured with respect to the interstellar medium 
in M87. Note also that the pressure in the knots might well exceed the ambient pressure 
without losing confinement, as long as the average pressure does not. 

We have calculated the average total pressure Pmean in the jet from the averaged 
emissivity and the best-fit 5-field for various angles and Lorentz factors, as shown in Figure 
^. The assumptions we have made in constructing this plot are analogous to those of |SB1V1| , 
who (arbitrarily) assumed a lower cutoff at 10''Hz, a high energy cutoff at lO^^Hz, a spectral 
index of a^o = 0.5, and equipartition between heavy-particle and electron energy. (Note: 
because the spectrum is steeper than = 0.5 above 10 GHz, our estimate of the pressure 
is likely to be an overestimate.) We have also assumed isotropic emission in the plasma rest 
frame by using an average value of ^9 = 54° for the term sin^^" ^ in the emissivity equation 
(0)- 

For comparison we have also calculated the equipartition pressure and plotted the 
ratio of pressure to equipartition pressure in Figure |^. It is obvious that we are far above 
equipartition for small values of Fjet and large 6'los, but as we approach the favored range 
of Fjet > 3 and 6'los ^ 30°, Ptotai approaches the equipartition value. The exact value of 
the pressure depends critically on the details we put into the model spectrum. For a jet 
composed entirely of electrons and positrons, the pressure would go down by a factor of 
i, whereas the equipartition pressure would only decrease by a factor of (|)7 = 0.67. The 
lack of information about the low-frequency spectrum inhibits any statements about the 
low-energy particle distribution. However, it is safe to assume that the power law does not 
continue down to non-relativistic energies. 
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It is obvious that a magnetic field far below equipartition alone cannot explain the 
behavior of the jet — it might account for the spectral changes but it requires the pressure 
to be much higher than that of the surrounding medium. Field orientations close to the 
line of sight will lead to an underestimate in emissivity, pressure, and intrinsic cutoff 
frequency and will require even lower magnetic fields and even higher pressures in order 
to prevent significant cooling. This changes as we increase Fjet: the inferred pressures are 
close to the value for the interstellar medium in M87 as derived by |White fc Sarazin (1988) , 
for Fjet ~ 3 — 5 and ^los ^ 25°. This, combined with the possible overpressure of the 
jet's immediate environment relative to the ISM, leads to the conclusion that there is no 
confinement problem. 

For small values of the magnetic field, one might ask if inverse Compton losses become 
dominant. A simple order of magnitude estimate shows that this is not the case. The ratio 
of synchrotron to inverse Compton loss timescales is equal to the ratio of photon energy 
density to magnetic field energy density ( [Rybicki fc Lightman 1979| ). The magnetic field 
strength has been estimated above. To derive an estimate of the photon energy density 
produced by the synchrotron emission, we normalize to the radio luminosity corrected for 
beaming and projection effects and integrate over an oro = 0.65 radio-to-optical power 
law that cuts off at 10^^ Hz. This shows that the inverse Compton lifetime due to just 
the synchrotron radiation field of the jet is roughly an order of magnitude longer than the 
synchrotron lifetime for the parameter range we suggested above. The starlight background 
at the center of M87 also contributes to the photon energy density. Using an isothermal 
sphere profile, normalized to the total luminosity of M87, we arrive at a central photon 
energy density roughly an order of magnitude smaller than that of the magnetic field, small 
enough to justify the assumption of neghgible Compton losses. 
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4.2. Comparison to Earlier Models 

It is instructive to compare our model to an earlier ad-hoc model by Meisenheimer 
et al. (1996, see also piretta fc Meisenheimer 1993| ), which bears a lot of similarity to 



our model. They start from the same assumption that the spectral changes along the 
jet can be explained by simple compressions and assume that the cutoff momentum 7c 
is almost constant along the jet, parameterizing it as a function only of the transverse 
jet radius (measured from the 2 cm radio map): 7c oc — note that for an adiabatic 

2 

compression of the plasma transverse to the flow 7c goes as R~i. They take the S-field 
to consist predominantly of a toroidal component, 5$, and hold the poloidal component 
Bz fixed. They determine the longitudinal compression ratio of -B$ and of the particle 
density n from their fit to the cutoff frequency v^. with equation p. However, in an 
adiabatic compression, the cutoff momentum varies as and will therefore be affected by 
longitudinal compressions as well (here is where our assumption of a disordered field allows 
us to determine a relation between density and magnetic field, so we can solve equation 
P uniquely for B). They neglect the fact that the synchrotron emissivity is enhanced in 
adiabatic compressions by n^+i" ■ (equation |T]) rather than n ■ Since the shocks 

might well be oblique, their assumptions that Bq> scales as the longitudinal compression 
ratio and that B^ is constant might also not be valid. 

Meisenheimer et al. (1996) favor an intrinsically onesided, subrelativistic jet, viewed 
close to perpendicular {Olos ~ 90°). Knot A would be a head-on shock in this scenario. 
As we have mentioned above, for this set of parameters additional acceleration has to be 
provided to maintain the optical emission out to large distances from the core. Meisenheimer 
et al. (1996) favor an unknown global acceleration process to explain the constancy of the 
cutoff momentum. With these assumptions, their model yields similar results to ours in 
that it reproduces the small scale brightness variations on the basis of the changes in cutoff 
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frequency. Our model could thus be regarded as an extension of their approach, putting it 
on the theoretical basis of adiabatic expansions, with a different mechanism for providing 
the large scale constancy of the spectrum. 

4.3. X-ray Emission 

An important result from the analysis above is that cooling longward of UV energies 
is not important over the length of the jet — the proper time is reduced by a significant 
factor and the magnetic fields are small enough to leave the spectral shape unchanged. This 
conclusion begs the question of the origin of the X-ray emission detected by the Einstein 
Observatory and ROSAT. Both observations show emission from the core/knot D region 
and knot A. 

4.3.1. The Role of Particle Advection 

Ultra-high-energy particles, capable of radiating in the X-ray regime, could be carried 
out from knot D, where X-ray emission is observed, to knot A, and reaccelerated in the 
shock by the adiabatic compression mechanism discussed above. For this to happen the 
spectrum would have to break rather than cut off in the optical. The presence of a break 
instead of a cutoff would not change our fits, as long as the break is located above the 
frequency we fitted. We have calculated the behavior of particles with X-ray emitting 
energies along the jet and found that for our best fit S-fields cooling out to knot A will 
have produced a spectral cutoff at ~ 10^^ Hz - which is where most of the Einstein HRI's 
sensitivity lies. Since the i?-fields we derived are upper limits, a lower field could leave 
the distribution function unchanged even at such high energies. Therefore this mechanism 
of producing X-rays is marginally consistent with our model. It might also account for 



at least part of the X-rays. Note that, as suggested by various authors (e.g., [Biretta ^ 
Meisenheimer 1993| ), the X-ray emission could also be of non-synchrotron origin altogether. 



4-3.2. Fermi Acceleration at Knot A 

The compression of a factor ~ 3 inferred from the analysis above indicates that Fermi 
acceleration might be present in knot A, although the moderate change in emissivity and 
the constancy in radio and optical spectral index suggests that it might not be very efficient. 
It is possible that Fermi acceleration occurs at parts of the shock only, resulting in a particle 
distribution superposed of a compressed and a Fermi-accelerated preshock distribution. 

In order for Fermi acceleration to take place at all the shock has to be subluminal, i.e., 
the intersection point of a given magnetic field line and the shock front has to move with 
a speed smaller than the speed of light, in which case we can find a frame in which the 
magnetic field is perpendicular to the shock front. This is the case for fields not too closely 
aligned with the shock plane. In the nonrelativistic case the field orientations leading to a 
superluminal shock are rare, and subluminal shocks are the rule rather than the exception, 
so one would expect Fermi acceleration to take place. 

Because in relativistic shocks the percentage of superluminal field orientations rises 
sharply with Fghock, Fermi acceleration should become less important. In this limit, most of 
the particle acceleration would occur through the mechanism of "shock drift acceleration" . 
|Begelman fc Kirk (1990)] presented a theory of this process valid in the relativistic case. 
They show that the adiabatic approximation is still accurate in the limit of Fp/Jp <^ 1, 
where Fp is the upstream Lorentz factor in the perpendicular shock frame and f3p the 
corresponding velocity. We have calculated this quantity for various obliquities and field 
orientations appropriate for knot A and it seems to fall into the desired range. For a 
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disordered magnetic field we would have to average the resulting spectrum over all possible 
field orientations. The more superluminal the shock the less important would effects of 
Fermi acceleration be. Relativistic corrections to shock drift acceleration are only important 
for high TpjSp, which, again, depends on the field orientation. Averaging over all possible 
field orientations would probably render these corrections unimportant. 

Fermi acceleration both changes the shape of an incoming power law spectrum and 
amplifies it. For an incoming electron spectrum of the form f{p) = AqP'"" in the test 
particle limit (i.e., the pressure provided by the accelerated particles is negligible), and 
a nonrelativistic shock, the change of the spectrum depends on s = where r is the 
compression ratio. If s < a the spectral index is changed to s. If s > a the slope remains 
unchanged but the spectrum is still amplified by a factor ( [Kirk 199^ . The radio 



spectral index is a ~ 0.5 — 0.65, implying a ~ 4 — 4.3. To provide a boost in emissivity by a 
factor of ~ 11, s has to be 5.1 — 5.5, implying a compression ratio of ~ 2.2 — 2.5 (assuming 
C = |). Note that this is very close to the compression ratio one derives for an adiabatic 
compression. The spectral index produced by such a shock is a ~ 1 — 1.25, consistent with 
the observed Oopticai ~ 1-2. Prury, Axford, fc Summers (1982)| showed that the produced 



power law softens as one departs from the test particle limit. Also, the simple treatment 
stated in this paragraph breaks down in the case of relativistic shocks, where the spectral 
index is no longer a simple function of the compression ratio. 



Recent investigations by [Ballard fc Heavens (1992")| have shown that oblique relativistic 



shocks can produce rather steep spectra, but other results indicate that they might be more 
efficient in accelerating particles (i.e., producing flatter spectra) than their nonrelativistic 
counterparts ( [Kirk fc Heavens 1989| ). [SBM[ find an optical-to-X-ray spectral index 
of aox ~ 1-4 for the knot A region, which seems to be consistent with low-efficiency 
acceleration. 
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If Fermi acceleration were present at the shock and effective enough to change the 
spectral shape in the optical, it would no longer be feasible to use the data from the whole 
jet to determine the magnetic field. Rather, the same analysis could simply be carried out 
seperately for the pre- and post-knot A regions of the jet. We would then have to make 
an estimate of the shock strength based on the known parameters in order to determine 
the ratio of pre- to post-shock Lorentz factor F. Even in this case we would need Lorentz 
factors of order Fjet ~ 3 to solve the cooling problem. 

However, based on the observed mild spectral changes, and the assumption that pitch 
angle scattering is strong, we conclude that Fermi acceleration, if present, will not be 
efficient enough to affect the spectrum below the cutoff. Inefficient Fermi acceleration might 
very well be present in knot A, producing the X-ray emitting particles observed. Prediction 
of the produced high-energy spectral index has to wait for more conclusive results on Fermi 
acceleration at relativistic oblique shocks. 



4.4. Polarization 

An important complication to the treatment above is the fact that the magnetic 
field will not be completely disorganized — polarization measurements show that in some 
regions a homogeneous component is present. In fact, it is possible that cancellations 
between regions with homogeneous fields but different orientations occur along the line of 
sight (pVIeisenheimer 1992|) . In such a case, the assumption of disorganized fields, leading to 
C |, is no longer justified. However, since the impact that the parameter ( has on the fit 
is minor (Fig. ^), we feel this caveat is not very severe and merely mention this complication 
here. 

In addition to the unknown orientation of the jet itself, the field orientation is also 
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unknown. Since synchrotron emission depends on the magnetic field orientation to the fine 
of sight 'd as sin^"*""" -(9, it can be strongly peaked away from the field direction. The cutoff 
frequency also depends on sind. Because the cutoff frequency is determined by sampling all 
regions along the line of sight, it is not obvious which value to choose for d. Fortunately, in 
a domain of disordered field the regions with the field oriented close to perpendicular will 
contribute most of the fiux, thus the error we make by setting sind to 1 will not be too 
large. Note that we used a value oi 'd = 53° in calculating the pressure, the appropriate 
average of sin"*^'^ over 47r sterad. 

Knot A shows a polarization of order 35% and a field orientation close to the shock 
plane. Neglecting relativistic effects, the small amount of upstream polarization, and the 
fact that shear might reduce the compression of the field, we can use the approximate 
formula ( |Hughes fc Miller 1995| ) 

^ a + \ 2-(l-r-2)cos2e' ^ ' 

where 11 is the fractional polarization, r is the proper compression ratio, and e is the line of 
sight angle from the shock plane, to obtain a lower limit on the compression ratio in knot A 
of r 2, valid for the inferred range of viewing angles with respect to the shock plane (BB). 
This is consistent with the compression ratios inferred above. Equation (0) also provides an 
upper limit of 35° on e; however, the shock is assumed to be oblique, so we cannot use this 
result to constrain Olqs- 



4.5. Particle Acceleration Radius 

Having an estimate of both magnetic field and cutoff momentum at the injection 
radius tq, we can now try to determine where the actual particle acceleration has to occur. 
We assume some radial dependence for the magnetic field in the inner portion of the jet 
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(i.e., smaller than 0.5"), for example a power law: B oc r~'^. Furthermore we make the 
simplifying assumption of a constant Fjet- By also taking B to be proportional to g'^, which 
determines the radial dependence of g, we can invert equation and solve for the radius 
at which the cutoff momentum approaches infinity, in other words, the minimum radius 
inside of which acceleration has to occur: 

race = ro(l + -r^) (8) 



A To 

where A = tr'^'^r 4 ; cind A = ,^ , A — -. The acceleration radius race approaches zero if 
A ■ ro — — 1. It is obvious that the estimate of raee depends critically on the value of a. In 
Figure || we have plotted race as a function of a for 6'los = 25° and Fjet = 3. In the same 
figure we have also plotted race for the case in which B no longer scales like g^ — in this 
case we have taken g oc r~^ and used the same values for 6'los and Fjet (dashed line). 

If the jet expands at constant opening angle and with uniform Fjet ? the decline of B 
with radius should correspond to a < 2, since the density scales like r~^ and dissipative 
effects will probably limit the rate of decline. Adopting an upper limit on the magnetic 
field at 0.01 pc of -B < 0.1 G ([Reynolds et al. 1996a|) limits cr to values smaller than 1. 



For Fjet = 3, cr = 1, and C = this constrains the acceleration radius to be race < 10 pc, 
or 0".06. On the other hand, the radius of acceleration cannot lie inside the Schwarzschild 
radius of the central black hole, which is of order 10~^ pc for a 10^ M© black hole. 

The dependence of race on a, (, and the core magnetic field is too strong to make 
any detailed predictions about where the acceleration actually has to occur. However, the 
race^curve is rather flat throughout most of the possible range for cr, which suggests that 
the most plausible value for faec falls between 1 and 10 pc. This is intriguingly far away 
from the central engine. 
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5. CONCLUSION 

We have proposed that the apparent lack of synchrotron coohng in the M87 jet 
most hkely indicates the presence of a sub-equipartition magnetic field. While the total 
(particle + magnetic) pressure needed to explain the observed synchrotron emissivity is 
uncomfortably high in the nonrelativistic limit, Doppler beaming effects consistent with 
bulk Lorentz factors in the range 2 — 5 lower the pressure requirements considerably. 
Fluctuations in the synchrotron emissivity and spectral cutoff frequency are consistent 
with adiabatic changes in the magnetic field strength and particle energies that accompany 
compressions and rarefactions along the flow — Fermi acceleration along the flow is not 
necessary to explain the observations. 

The knots are identified with relatively weak shocks, as inferred from other data by 



BB| . The first-order Fermi acceleration expected to occur at such shocks, if any, would 
generate a particle energy distribution steeper than the n{E) oc E^'^ needed to produce the 
radio-to-optical synchrotron spectrum. Thus, effects of particle acceleration along the jet 
might be apparent only shortward of the cutoff frequency, e.g., in the X-ray band. However, 
a disordered magnetic field or a "superluminal" field orientation with respect to the shock 
front (Begelman & Kirk 1990) could further reduce the efficiency of Fermi acceleration, 
hence we should continue to regard the origin of the X-ray emission as unknown. 

The pressure estimates we derive are consistent with the assumption that the M87 jet 



is embedded in a moderately overpressured bubble, as suggested by |BB| . As a result, it 
seems that the set of parameters we have suggested above can solve the cooling problem of 
the M87 jet, as sub-equipartition fields are able to explain both the behavior of the cutoff 
frequency and the confinement of the jet. 

Using derived values for the magnetic field and the cutoff momentum at ro we can put 
an upper limit of 10 pc on the radius at which most of the particle acceleration occurs. 
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Explaining why the acceleration is confined to a particular scale (which may be quite large 
compared to the size of the central black hole) poses an interesting problem for future work. 

The method we have developed should be apphcable to other radio sources as well, 
as long as propagation effects on the particle momentum distribution can be treated as 
adiabatic. Thus, as soon as comparable data becomes available for other jets, it would be 
easy to apply the same analytic techniques to them and obtain limits on the magnetic field 
and the pressure in these sources. 

The analysis of the M87 jet we carried out in this paper could be improved by better 
deprojection models of the jet. Higher quality spectral data both in the near infrared and 
the X-ray would help to verify the assumption of constancy of the underlying spectral 
shape and decrease the uncertainty in the location of spectral features (the knowledge of 
which is essential for this technique to work). HST/NICMOS will be ideally suited to 
mapping out the jet properties in the near infrared, a spectral region in which the run of the 
spectral index is still fairly poorly determined. AXAF, providing imaging spectroscopy at 
sub-arcsecond resolution, will be the ideal tool to obtain more detailed information about 
the X-ray spectrum and emissivity along the jet, hopefully allowing us to determine the 
radiation mechanism at high energies. 

The model we have described above fits in nicely with the new picture of M87 
that has emerged in several recent papers. A moderately relativistic jet and magnetic 
fields somewhat below equipartition can explain the lack of cooling observed at optical 
wavelengths. Oblique shocks, weak enough to render Fermi acceleration unimportant, can 
explain both the small scale variability of the cutoff frequency and the brightness variations 
via adiabatic compressions. We believe this model can tie together most of the observations 
available today and will take us closer to the true nature of the jet. 
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Fig. 1. — Measured cutoff frequency along tlie jet from NMR| (vertical lines with error bars). 



Solid line: best fit curve as calculated from the emissivity measured by ^BM| for Fjet = 1-1, 
C = |, and 6'los = 90°; dashed line: best fit curve for the case of equipartition calculated 



from 3BM data for the same set of parameters. 



Fig. 2. — Contour plot of the values as a function of -Bmcan and the injection cutoff 
momentum Pq for Q = 25° and Fjct = 3.0. Shown are the contours corresponding to a 
reduced of 45 and integer multiples of this value. The shaded line on the right of the plot 
shows the equipartition S-field for this set of parameters. 

Fig. 3. — 5-field averaged along the jet between 0".5 and 22" as a function of Fjet for 
^Los = 15° (solid line), ^los = 20° (dotted line), ^los = 25° (dashed line), and ^los = 30° 
(dashed-dotted). The grey region indicates the jet-to-counterjet brightness ratio limit of 
150 (dashed grey boundary) - 380 (solid grey boundary, piretta 1993|) . 



Fig. 4. — Ratio of the best-fit 5-field to the respective equipartition S-field B^q averaged 
between 0".5 and 22" for ^ = | as a function of Fjet- The different values of Olos are labeled 
similarly to Fig. (|^). The grey region shows the jet-to-counterjet brightness limit (dashed 
grey boundary: 150, solid grey boundary: 380). 

Fig. 5. — Fractional deviation ^ of the best fit -B-field (averaged between 0".5 and 22" for 
a. two different values of C- = where Ci = | and either C2 = ^ (hatched, 

solid boundary) or ^2 = 1 (hatched, dotted boundary); b. a uniform jet, compared to a 
jet slowing down at knot A: ( ^) = ^'^"'"0"""'^'"'°''''^ , where -Bbrcak is the best fit average 

\ / p ^uniform 

S-field for a jet slowing down from Tj^t^A- to Tjet,A+ at knot A, and -Buniform is the best fit 
average field for a uniform jet with Fjet = ^jet,A- (dark grey region, short dashed boundary) 
and Fjet = ^jet,A+ (light grey region, long dashed boundary). The latter is plotted versus 
^jet,A+'i c. two different 2cm radio spectral indices: (■^) = 1^^^ b^^I standard 
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value «! = 0.5 and a2 = 0.65, shown as the black area. (The width of the band corresponds 
to ^Los between 15° and 30° in each case.) 

Fig. 6. — Total pressure Pmean averaged along the jet between 0".5 and 22" as a function of 
Fjet. Labels according to Fig. (0). For comparison, the dashed-triple-dotted line shows the 
equipartition pressure for a non-relativistic jet seen edge on (i.e., 6'los = 90°). The hatched 
area shows the estimated ISM pressure (White & Sarazin 1988) in M87. 

Fig. 7. — The ratio of the best-fit particle pressure to the respective equipartition value 
averaged along the jet as a function of Fjet. Labels according to Fig. (H). 

Fig. 8. — The minimum acceleration radius race, as a function of a (5 oc r'^) for Fjet = 3 
and 6'los = 25° in the case oi B oc g'^ (solid line). The dashed line shows race for the case of 
g oc r~^ and the same values of Fjet and 6'los- The hatched area indicates the limit set by 
10 Schwarzschild radii for the ~ 10^ Mq central black hole. The dashed-dotted line shows 
the (arbitrary) injection radius tq = 0".5 or 80 pc. 
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